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Abstract 



The interaction potential V(Q ) between static test charges can be used 



^O ', to define an effective charge ay(Q ) and a physically-based renormalization 



scheme for quantum chromodynamics and other gauge theories. In this paper 

we use recent results for the finite-mass fermionic corrections to the heavy- 

c | ■ quark potential at two-loops to derive the next-to-leading order term for the 

&h, Gell Mann-Low function of the V-scheme. The resulting effective number of 

Q_i! flavors Np(Q 2 /m 2 ) in the ay scheme is determined as a gauge-independent 

and analytic function of the ratio of the momentum transfer to the quark 

pole mass. The results give automatic decoupling of heavy quarks and are 

independent of the renormalization procedure. Commensurate scale relations 

then provide the next-to-leading order connection between all perturbatively 

calculable observables to the analytic and gauge-invariant ay scheme without 

any scale ambiguity and a well defined number of active flavors. The inclusion 

of the finite quark mass effects in the running of the coupling is compared with 

the standard treatment of finite quark mass effects in the MS scheme. 
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I. INTRODUCTION 

Is there a preferred effective charge which should be used to characterize the coupling 
strength in QCD? In principle, any perturbatively calculable observable can be used to define 
an effective charge. In quantum electrodynamics, the Dyson running coupling a QED (Q 2 ), 
defined from the linearized potential between two infinitely-heavy test charges, has been used 
as the traditional coupling. The corresponding definition of the non-abelian QCD coupling 
is customarily given by identifying the ground state energy of the vacuum expectation value 
of the Wilson loop as the potential V between a static quark- ant iquark pair in a color singlet 
state 0: 

V(r,m 2 ) = - hm - log (0 1 Tr (p exp U dx^T 11 ) j |0) (I) 

where r denotes the relative distance between the heavy quarks, m is the mass of "light" 
quarks contributing through loop effects, and T a are the generators of the gauge group. It 
is then convenient to define the effective charge ay(Q 2 , m 2 ) as 

K(QW) ^ -4nC F ^f^ (2) 

in momentum space. The factor Cp is the value of the Casimir operator T a T a of the 
external sources (which are in the fundamental representation) and factors out to all orders 
in perturbation theory, and Q 2 = —q 2 is the spacelike momentum transfer. 

The effective charge ay{Q) provides a physically-based alternative to the usual modified 
minimal subtraction (MS) scheme. As in the corresponding case of Abelian QED, the scale 
Q of the coupling ay(Q) is given by the exchanged momentum. There is thus no ambiguity 
in the interpretation of the scale. All virtual corrections due to massive fermion pairs are 
incorporated in ay through loop diagrams which depend on the physical mass thresholds. 
When continued to time-like momenta, the coupling has the correct analytic dependence 
dictated by the production thresholds in the crossed channel. Since ay incorporates quark 
mass effects exactly, it avoids the problem of explicitly computing and resumming quark 
mass corrections which are related to the running of the coupling. Thus the effective number 
of flavors Np(Q/m) is an analytic function of the scale Q and the quark masses m. The 
effects of finite quark mass corrections on the running of the strong coupling were first 
considered by De Riijula and Georgi [fj] within the momentum subtraction schemes (MOM) 
(see also [§~|§]). The two-loop calculation was first done by Yoshino and Hagiwara j?| in the 
MOM-scheme using Landau gauge and also recently by Jegerlehner and Tarasov || using 
background field gauge. 

One important advantage of the physical charge approach is its inherent gauge invari- 
ance to all orders in perturbation theory. This feature is not manifest in massive /3-functions 
defined in non-physical schemes such as the MOM schemes. A second, more practical, ad- 
vantage is the automatic decoupling of heavy quarks according to the Appelquist-Carazzone 
theorem [HI. 



By employing the commensurate scale relations |L0| other physical observables can be 



expressed in terms of the analytic coupling ay without scale or scheme ambiguity. The 
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quark mass threshold effects in the running of the coupling are taken into account by utiliz- 
ing the mass dependence of the physical ay scheme. In effect, quark thresholds are treated 
analytically to all orders in m 2 /Q 2 ; i.e., the evolution of the physical ay coupling in the in- 
termediate regions reflects the actual mass dependence of a physical effective charge and the 
analytic properties of particle production. Furthermore, the definiteness of the dependence 
in the quark masses automatically constrains the scale Q in the argument of the coupling. 
There is thus no scale ambiguity in perturbative expansions in ay. 

In the conventional MS scheme, the coupling is independent of the quark masses since the 
quarks are treated as either massless or infinitely heavy with respect to the running of the 
coupling. Thus one formulates different effective theories depending on the effective number 
of quarks which is governed by the scale Q; the massless /^-function is used to describe the 
running in between the flavor thresholds. These different theories are then matched to each 
other by imposing matching conditions at the scale (normally the quark masses) where the 
effective number of flavors is changed. The dependence on the matching scale can be made 
arbitrarily small by calculating the matching conditions to high enough order. For physical 
observables one can then include the effects of finite quark masses by making a higher-twist 
expansion in m 2 /Q 2 and Q 2 /m 2 for light and heavy quarks, respectively. These higher-twist 
contributions have to be calculated for each observable separately, so that in principle one 
requires an all-order resummation of the mass corrections to the effective Lagrangian to give 
correct results. 

The specification of the coupling and renormalization scheme also depends on the defi- 
nition of the quark mass. In contrast to QED where the on-shell mass provides a natural 
definition of lepton masses, an on-shell definition for quark masses is complicated by the 
confinement property of QCD. In this paper we will use the pole mass m which has the 
advantage of being scheme and renormalization-scale invariant. 

A technical complication of massive schemes is that one cannot easily obtain analytic 
solutions of renormalization group equations to the massive (3 function, and the Gell-Mann 
Low function is scheme-dependent even at lowest order. 

In this paper we present a two-loop analytic extension of the ay-scheme based on the 



recent results of Ref. [TT|. The mass effects are in principle treated exactly to two- loop 
order and are only limited in practice by the uncertainties from numerical integration. The 
desired features of gauge invariance and decoupling are manifest in the form of the two-loop 
Gell-Mann Low function, and we give a simple fitting-function which interpolates smoothly 
the exact two-loop results obtained by using the adoptive Monte Carlo integrator VEGAS 



[ 12 1 . Strong consistency checks of the results are performed by comparing the Abelian limit 
to the well known QED results in the on-shell scheme. In addition, the massless as well as 
the decoupling limit are reproduced exactly, and the two-loop Gell-Mann Low function is 
shown to be renormalization scale (/i) independent. 

As an application we show how the analytic ay-scheme can be used to calculate the 
non-singlet hadronic width of the Z-boson, including finite quark mass corrections from the 
running of the coupling and compare with the standard treatment in the MS scheme where 
the corresponding effects are calculated as higher twist corrections. 

Recently (see Ref. [0) we proposed an alternative way of incorporating mass effects 
connected with the running of the coupling by making an analytic extension of the MS- 
scheme where the coupling is an analytic function of both the scale and the quark masses. 



This analytic extension of the MS-scheme is defined by connecting the MS coupling to the 
^-scheme using a commensurate scale relation based on the (BLM) scale-setting procedure 



14| . The new modified coupling a-^iQ) inherits most of the good properties of the ay 
scheme, including its correct analytic properties as a function of the quark masses and its 
unambiguous scale fixing |~3] . 



However, the conformal coefficients in the commensurate scale relation between the ay 
and MS schemes do not preserve one of the defining criterion of the potential expressed in 
the bare charge, namely the non-occurrence of color factors corresponding to an iteration 
of the potential. This is probably an effect of the breaking of conformal invariance by the 
MS scheme. The breaking of conformal symmetry has also been observed when dimensional 



regularization is used as a factorization scheme in both exclusive [|T|,|nj and inclusive [17 
reactions. Thus, it does not turn out to be possible to extend the modified scheme 5m§- be- 
yond leading order without running into an intrinsic contradiction with conformal symmetry. 
Note, however, that this difficulty does not affect using the MS scheme as an intermediate 
renormalization scheme when connecting physical observables. For completeness we give the 
results of such an extension in an appendix. 

The paper is organized as follows: In section [II] we derive the second term of the Gell- 
Mann Low function in the physical V^-scheme as a renormalization-scale-independent func- 



tion of the ratio of the physical momentum transfer Q and the pole mass m. In section [LI I 
we present numerical results of the effective number of flavors and compare it with results 
obtained in the gauge- dependent momentum subtraction schemes. In addition, various con- 
sistency checks are performed, and numerical fits are presented. In section [TV] we illustrate 
some of the properties of the analytic ay scheme and demonstrate the effect of the quark 
mass thresholds on the mass-dependent evolution and compare with the massless evolution. 
In section [V] we compare the calculation of the hadronic width of the Z-boson in the analytic 
ay scheme to the conventional MS scheme with a mass-independent coupling and explicit 
higher-twist corrections for mass effects. In section [VT] we summarize our results and indicate 
future applications. The definition of the analytically extended scheme <5ms beyond leading 
order is discussed in the appendix. 

II. THE GELL-MANN LOW FUNCTION THROUGH TWO LOOPS 

The physical charge ay(Q,m) can be expressed as a perturbative series in any other 
renormalization scheme. For example, in the minimal subtraction scheme, the perturbative 
series has the form: 

av(g,m) = a M8 (/x)fl + t; 1 (g,m(/x)^)^^ + ^(g,m(^),/i)^M + ..) (3) 

Y IT 7T Z J 

where the massless limit of the coefficients V\ and v 2 are known in the literature [fT|JT~8| |23 



Since the physical charge ay(Q,m) cannot depend on the renormalization scale /x, the in- 
dependence on the right-hand side of Eq. (fj) must cancel to the order we are working. 
Notice that the coefficients also depend on the renormalization scale /x used for the mass 
renormalization, i.e. through the dependence of the running mass m(/x). Fig. p] shows the 
Feynman diagrams for the fermionic contributions to the two-loop coefficient V2(Q, m(/x), /x). 



These contributions depend on the mass renormalization used for the one-loop coefficient 
Vi(Q , m(fi) , fi) . Since we are predominantly interested in the flavor-threshold dependence of 
heavy quarks, we shall relate the running mass to the pole mass which is renormalizat ion- 
scale independent and gives explicit decoupling. This also provides a physical picture as 
well as a straightforward Abelian limit. 

The next-to-leading order relation between the MS mass m(fi) and the pole mass m is 
given by |E3|, 



m(fi) = m 



1 - C F 1 + - log 7 - log 4vr 

7i V 2 m 4 



(4) 



where 7 is the Euler constant. Inserting Eq. (|]) into Eq. (||) gives at next-to-next-to-leading 
order 



a v (Q,m) = a MS (n) 



l + Vl (Q,m, / i)^^ + h(Q,m, / i) + A m (g,m^)] ^ (//i 
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7T 



(5) 



where A m (Q,m, fi) denotes the contribution arising from v\ when changing from the MS 
mass to the pole mass: V\(Q , m(fj.) , /x) = t>i(Q,m,/i) + A m (Q, m, /x) " M s w . 

The Gell-Mann Low function |25j for the V^-scheme is defined as the total logarithmic 
derivative of the effective charge with respect to the physical momentum transfer scale Q: 



qiy[ 9.) = da i(Q^ =f;-vff 



?» 



dlogQ 



i=0 



4 +2 (Q,m) 

7T'' 



i+1 



(6) 



where in the massless case the coefficients ipy and ipy are given by, 



,(0) 



^ ( m = 0) 



— N c N F = N F 

6 3 2 3 



^W ( m = 0) = — M - —N C N F _ -C F N F = — - —N F . 
rv v v 12 c 12 4 4 12 



(7) 
(8) 



For the massive case all the mass effects will be collected into a mass- dependent function 
Np. In other words we will write 



# i 9 - 



Q 



2 3 Fy \m 
51 19^(1) (Q 

T 12 ^ m 



(9) 
(10) 



where the subscript V indicates the scheme dependence of N FV and N FV . 

Taking the derivative of Eq. (|5|) with respect to logQ and re-expanding the result in 
av(Q,m) gives the following equations for the first two coefficients of \lV : 



4P (2 



dv l (Q } m,ij,) 

dlogQ 
d[v 2 (Q, m, fi) + A m (Q, m, //)] 
dlogQ 



+ 2fi(Q,m,/x) 



dvi(Q,m,n) 
dlogQ 



(11) 
(12) 
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FIG. 1. The two-loop massive fermionic corrections to the heavy quark potential (from |TT|), 
Double lines denote the heavy quarks, single lines the "light" quarks with mass m. The first two 
rows contain diagrams with a typical non-Abelian topology. The middle line includes the infrared 
divergent "Abelian" Feynman diagrams. They contribute to the potential only in the non-Abelian 
theory due to color factors ex CfCa- In addition, although each diagram is infrared divergent, 
their sum is infrared finite. The infrared finite Feynman diagrams with an Abelian topology plus 
the diagrams consisting of one-loop insertions with non-Abelian terms are shown in the last two 
rows. 



The argument Q/m indicates that there is no renormalization-scale dependence in Eqs. ( |TTD 
and Cp72|) . Rather, ipy and ipy are functions of the ratio of the physical momentum transfer 



Q = V~Q 2 an d the pole mass m only. The expression for ip v agrees with our result in 
Ref. JT3J. In Eq. ( |T2] ) the derivative of the A m (Q, m, /x)-term comes from using the pole-mass 



instead of the MS mass, whereas the remaining mass dependence in Eq. (|12]) is arbitrary in 
the sense that a different mass scheme is formally of higher order. In addition we note that 
the contribution 2vidvi/dlogQ cancels the reducible contribution (labeled 2vp in Fig. [1]) 
to t>2; it is thus sufficient to consider one quark flavor at a time. 

III. NUMERICAL RESULTS FOR THE ANALYTIC N F 

Because of the complexity of the integrals encountered in the evaluation |TT1 of the mas- 



sive two-loop corrections to the heavy quark potential, the results were obtained numerically 
using the adoptive Monte Carlo integration program VEGAS [12]. Thus the derivative of 
the two-loop term v% was calculated numerically, whereas the other terms in Eqs. ( [TT]) and 
([12]) were obtained analytically. The results are given in terms of the contribution to the ef- 
fective number of flavors N F y (-*A and N FV (Q) in the V-scheme, from a given quark with 
mass m, defined according to Eqs. (||) and ( |10|) respectively. The Appelquist-Carazzone [f| 
theorem requires the decoupling of heavy masses at small momentum transfer for physical 
observables. Thus we expect N FV (Q) to go to zero for Q/m — * 0. The massless result 

N FV — > 1 must also be recovered for large scales. 



The calculation presented in Ref. |TT| required the evaluation of four-dimensional inte- 
grals over Feynman parameters. Our results are based on 50 iterations of the integration 
grid each comprising 10 7 evaluations of the function which were needed to achieve adequate 
convergence. Even so, the Monte Carlo results still are not completely stable for small values 
of Q/m, especially in the light of the numerical differentiation required in Eq. ([12]) . Never- 
theless, accurate results can be obtained by fitting the numerical calculation to a suitable 
analytic function as shown below. 

The one-loop contribution to the effective number of flavors Np follows from the stan- 



dard formula for QED vacuum polarization. In our earlier paper [fL3|] we used the simple 
representation in terms of a rational polynomial |2|]: 

<> (9.)„-L^ = °- 2 g ,13) 

F > v \m) 1 + 5^ 1 + 0.2^ V ; 

which displays decoupling for small scales and the correct massless limit at large scales. 
Similarly, the numerical results for the two-loop contribution can be fit to the form 

ai 91 + a2 9l 

N (i) (Q\ l m? 2 m 4 ha\ 

^ Km) , , Q 2 , Q 4 { } 

m z m 4 

The parameter values a, and the errors obtained from the fit to the numerical calculation 
in the ^-scheme for QCD and QED are given in Table |. Similar decoupling forms have 
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TABLE I. Values for the parameters a\, 0,2 and 03 in the ^/-scheme for QCD and QED obtained 
by fitting our numerical results to the form given by Eq. (|l4|). The \ 2 values were obtained by 
subscribing a constant 0.01 error to each data point in the fit. (Here d.o.f. denotes the effective 
number of degrees of freedom for the fit, i.e. the number of fitted points minus the number of 
parameters). 



a\ a 2 03 

QCD -0.571 ±0.034 0.221 ± 0.015 1.326 ±0.116 



j£ 



d.o.f. 
86 
43 

QED 1.069 ±0.0088 0.0133 ± 0.0002 0.402 ± 0.005 ^ 



been used for interpolating the flavor dependence of the effective coupling in the momentum 
subtraction schemes (MOM) |7||| . 

In the case of QCD we obtain the following approximate form for the effective number 
of flavors for a given quark with mass m: 

(-0.571 + 0.221%)% 

AT« (9.\ « A m > m (15) 

Fy \m Q 2 Q 4 (15j 

1 + 1.326-+T + 0.221- 



m 2 m 4 



and for QED 



(1.069 + 0.0133%)% 

Fy \m O 2 O 4 ' [ } 

V y 1 + 0.402^- + 0.0133^- 



m 2 m 4 



The results of our numerical calculation of Np V in the ^/-scheme for QCD and QED 
are shown in Fig. || The decoupling of heavy quarks becomes manifest at small Q/m, and 
the massless limit is attained for large Q/m. The QCD form actually becomes negative at 
moderate values of Q/m, a novel feature of the anti-screening non-Abelian contributions. 
This property is also present in the (gauge dependent) MOM results. In contrast, in Abelian 
QED the two-loop contribution to the effective number of flavors becomes larger than 1 at 
intermediate values of Q/m. We also display the one-loop contribution N FV (Q) which 
monotonically interpolates between the decoupling and massless limits. The solid curves 
displayed in Fig. |2| shows that the parameterizations of Eq. ( |T5D which we used for fitting 
the numerical results are quite accurate. This is also indicated by the x 2 values obtained 
for the fits as given in Table [j| 

A strong check of our results, as well as the results presented in Ref . [O] , is the agree- 



ment with the two- loop Gell-Mann Low function for QED p6|-f2q1. Fig. |3] contains detailed 
comparisons of the analytic QED result with our numerical computation. For comparison 
the figure also displays the purely non-Abelian part of the QCD result as well as the total 
QCD result. The scalar functions occurring in the Abelian corrections are also used in the 
evaluation of the non-Abelian contributions (oc Ca), and it is therefore important to know 
that they were calculated correctly. 




FIG. 2. The numerical results for the gauge-invariant Npy in QED (open circles) and QCD 
(triangles) with the best x 2 fits of Eqs. ( |16D and ( |T5| ) superimposed respectively. The dashed line 
shows the one-loop N FV function of Eq. (|17|). For comparison we also show the gauge dependent 
two-loop result obtained in MOM schemes (dash-dot) |7]||. At large ^ the theory becomes effec- 
tively massless, and both schemes agree as expected. The figure also illustrates the decoupling of 
heavy quarks at small — . 



N 



(i) 2 

F,V 



1.5 



0.5 - 



-0.5 



- 








o 


:QED(a v ) 






A 


:QCD(a v ) 






A 


: QCD (CA) 1 






- 




.'' y 


' i A a ^ A 


- 




A 








1 





10 



10 



Q/m 



(i) 



FIG. 3. Comparison of the Abelian limit of our results (open circles) for Npy based on the 



calculation in Ref. 1 11] which was done in the MS-scheme with the well known result in the literature 



1 2(3 2£ ] done in the on-shell renormalization scheme (solid line) . Also shown are the gauge invariant 
non- Abelian contribution only (oc Ca) (open triangles) as well as the sum of all terms in QCD (solid 
triangles). The correct Abelian behavior is a very strong check on the results given in Ref. |lH| . All 
Monte Carlo results are based on 10 evaluations per iteration and 50 iterations of the integration 
grid. 
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FIG. 4. Illustration of the renormalization scale independence of the two-loop effective number 
of flavors N F v as a function of the ratio of the physical momentum transfer Q over the pole mass m. 
Numerical instabilities are visible for small values of ^ and occur because of limited Monte Carlo 
statistics (10 7 evaluations for each of the 50 iterations). The two fits obtained, which agree within 
statistical errors, are shown as a solid and dashed line for fi = m and fi = 0.031m respectively. 

Another important test of our results is renormalization-scale (/i) independence, which 
follows from the fact that the effective number of flavors in the ^-scheme is a physical 
quantity. This is illustrated in Fig. ^ which shows the results obtained for two different 
renormalization scales (/i = 0.031m and fi = m). The figure also shows the fits obtained for 
the two different cases. In fact, the differences are so small that the two lines cannot easily 
be distinguished. 

We can also apply the same fitting procedure to the dependence of the one-loop effective 
N F : 



N 



(o) 

F,V 



1 



X* 



1 + (5.19 ± 0.03) §r d.o.f. 



19 

27 



(17) 



This gives a higher precision global fit compared to the form in Eq. fll3|) . 



IV. SOME PROPERTIES OF THE ANALYTIC COUPLING IN THE V-SCHEME 

Using the numerical results for N FV and N FV the evolution equation (|[) can be 
solved numerically using the classical Runge-Kutta algorithm. As starting value we use 
av{Mz,m) = 0.126 in next-to-leading order and av{Mz,m) = 0.134 in leading order which 
have been obtained from the value aig(Mx) = 0.118. It should be noted that it is straight 
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forward to solve this equation numerically since we are using the pole-masses which do not 
depend on Q. This should be compared with the MOM-scheme where one gets two coupled 
differential equations to solve, both for the coupling and the mass. 

The resulting leading and next-to-leading order \l/-function in the V-scheme is shown 
in Fig. ^| scaled with the leading dependence on ay, i.e. — ^ v I {oty / tt) . For comparison 
the figure also shows the ^/-function obtained with discrete theta-function thresholds with 
continuous matching, av(Q,Np = 0) = av{Q,Np = 9 + 1), at the naive matching scale 
Q = m. As can be seen from the figure there are significant differences between the two 
approaches both in leading and next-to-leading order. In fact the difference becomes larger 
when going to next-to-leading order. We also note that the scale dependence of the coupling 
is larger in next-to-leading order and that the convergence of the \l/-function is not very 
good for scales below a few GeV. From the figure it is also clear that there are no plateaus 
in the analytic treatment of quark masses. Thus there is no region of the scale Q below ~ 1 
TeV where all quark masses can be neglected at the same time. 



-y v /(a2/7t) 




10 10 

Q [GeV] 



FIG. 5. The scaled ^-function, -$v , /(tty/"') in the analytic V-scheme av(Q,rm) (solid) com- 
pared to the ay(Q,0) scheme with discrete theta-function treatment of flavor thresholds with 
continuous matching at Q = m (dashed). 

The solution of the evolution equation also gives the coupling as a function of the scale 
Q. The relative difference between the analytic av(Q,mi) and the discrete theta-function 
treatment of flavor thresholds with continuous matching at Q = m, ay{Q, 6), is shown in 
Fig. ^ both in leading and next-to-leading order. 

As can be seen from the figure, the difference between the analytic and step-function 
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[a v (Q,mj) - a v (Q,0)]/a v (Q,m i ) 
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Q [GeV] 

FIG. 6. The relative difference between the solutions to the evolution equation using the 
analytic ^-function, av(Q,rrii), versus discrete theta-function thresholds, ay(Q,G). The solid 
(dashed) curve shows the next-to- leading (leading) order result. 

treatment of quark masses in the running persists when going to higher order. In fact we 
expect this difference to remain to all orders. The reason is that the \1/ function is not 
continuous in the step-function approach and the stepsize at the thresholds is governed 
by the lowest order term ip(°\ Thus there will always be a finite difference between the 
continuous \l/-function and the one with theta-function thresholds. The difference can be 
made smaller by modifying the matching scale to be Q = 3m (but still using continuous 
matching) which is also illustrated in the figure. However, the difference cannot be made 
smaller than ~ 1%. The only way to include the finite quark mass effects in the fixed flavor 
treatment is by making a higher twist analysis to all orders in m 2 /Q 2 and Q 2 /m 2 for light 
and heavy quarks respectively. 

Noting that the differential equation for the scale dependence of the coupling is homo- 
geneous in Q/m we can also get the logarithmic derivative of ay(Q, rrii) with respect to the 
heavy quark masses, 



dav(Q,rrii 
dlogm 



dav(Q,rrii 



dlogQ 



(18) 



where the subscript q refers to the quark part of the \l/-function. The resulting mass- 
dependence is shown in Fig. [I] for each of the heavy quarks. 

The figure illustrates how the quarks decouple for small scales Q and how they become 
effectively massless for large scales Q. In the intermediate region the mass- dependence 
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(d oc v / d log m) / (1/3 all% + 19/12 al/ic) 
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FIG. 7. The logarithmic derivative of ay(Q, mi) with respect to the quark mass for each of the 
heavy quarks. 

depends on the ratio of the mass to the scale Q. For large scales Q the derivative approaches 
the asymptotic value 



day (Q, rrii 



dlogm 



L ~S>m Jj 



la v (Q,mi) 5Nc + 3C F a v (Q,mi 
3 



n 



12 



7T 



(19) 



which is the same as the Np dependent part of the \l/-function apart from the differing sign. 
This can also be derived from the decoupling relations for matching fixed Np couplings as 
for example is done in the MS scheme. 

V. APPLICATION 



The purpose of this section is to compare the treatment of finite quark mass effects in 
the V-scheme with the standard treatment in the MS scheme. To do this comparison we 
will follow our earlier paper |L3| and use the non-singlet hadronic width of a Z-boson with 
arbitrary mass yfs starting from the physical mass ^/s = Mz for normalization. 

The finite quark mass effects that we are interested in are in leading order given by the 
"double bubble" diagrams, which are shown in Fig. |], where the outer quark loop which 
couples to the weak current is considered massless and the inner quark loop is massive. These 
corrections have been calculated in the MS scheme as expansions in m? / s [fj9| and s/ttiq |3^ 



for light and heavy quarks, respectively, whereas they have been calculated numerically in 
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1 31 1 . In. addition the a^ correction due to heavy quarks has been calculated as an expansion 



in s/toq in [32]. Other types of mass corrections, such as the double-triangle graphs where 



the external current is electroweak, are not taken into account. 






FIG. 8. The "double bubble" diagrams. The crosses represent the external electro weak current, 
the thin line is a massless quark and the thick line is a massive quark. 



The non-singlet hadronic width of a Z-boson with arbitrary mass y/s is given by 



had\ 



^Vew) 2 +(^) 2 } 



1 + i Cf ^t- 



(20) 



where a^ s (s) is the effective charge |33| which contains all QCD corrections. In the following, 
the next-to-leading order expressions for the effective charge a^ s (s) in the MS and V schemes 
will be compared for arbitrary s using next-to-leading order evolution starting from the 
physical mass y^ = Mz which is used as normalization condition. 



A. MS scheme treatment 



In the MS scheme the effective charge a^ s (s) is given by 
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where the coefficients T\ and r2 are given by, 

ri,Ms(/i = y/s) = --C F + — N c + (— - 2C 3 ) A, 

= 1.986 -0.115iV F 
r-N 23^ 2 / 53 1L\ r2 / 101 11. , ., ,.. 

^<" = ^ = "£<* + ("144 " T C ') N ° + (-192 + T & I CfJV ° + 
11 >\ „ /151 19 > tt 2 \ „„ 



(21) 
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(with /3 = Vy ( m = 0) anc ^ A = W ( m = 0)) and the functions F and G are the effects of 
non-zero quark masses for light and heavy quarks, respectively. The expansions of the a 2 ^ 
finite quark mass corrections are given by[] 
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(22) 
(23) 



which are accurate to within a few percent for m 2 / s < 0.25 and s/m 2 n < 4 respectively. We 



will also use the relation |3l 

,2 





11 2 A 

12 3 s 



(24) 



to obtain F\ in the interval 0.25 < m 2 /s < 1 where the expansion of F\ given above breaks 
down. 

The a|jg finite quark mass corrections are only known for heavy quarks (G2), whereas 
the corresponding corrections due to light quarks (F2) have not yet been calculated. The 
known corrections in G2 are small, for itiq = s they are of order G<i ~ 0.1 for the top quark. 

The number of light flavors N L in the MS scheme is a function of the renormalization 
scale /1. In the following we will assume that the matching of the different effective theories 
with different number of massless quarks is done at the quark masses. In other words a quark 
with mass m < \x is considered as light whereas a quark with mass m > // is considered as 
heavy. In addition the MS quark masses are used. The general matching condition in the 
MS scheme is to next-to-leading order given by ||34|| , 
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(25) 



where m(fi) is the mass of quark number Nl + 1. The dependence on the matching scale 
can be made arbitrarily small by calculating the matching condition to high enough order. 
However this does not mean that the finite quark mass effects are taken into account. The 
only way to include these mass effects in the ordinary MS treatment is by making a higher 



x In our earlier paper [13] there was a typographical error giving the wrong signs for the two 
ln-terms in G\. 



15 



twist expansion to all orders in m 2 /Q 2 and Q 2 /m 2 for light and heavy quarks respectively, 
i.e. the functions F and G given above. 

In the following comparison we will restrict ourselves to the next-to-leading order ex- 
pression for a^ s {s) in the MS scheme including the finite quark mass corrections, i.e. 



a#°(s) = °&W + 



MS 
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with // = -y/i and next-to-leading order matching done at the quark masses. 



B. V scheme treatment 



(26) 



In order to relate the hadronic width of the Z-boson to the ay scheme we will require 
the massless coefficients in the relation between a^§ and ay for Q 2 ^> m 2 , 
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where the coefficients are given by 
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The logarithmic Qj \x dependence of the coefficients follows from requiring the expansion of 
ay(Q) in a-^(fi) to be /i independent to the order we are working. Inserting Eq. fl27D into 
the massless version of Eq. ( pT|) (i.e. without the finite quark mass corrections F and G 
for light and heavy quarks respectively) gives the relation between the effective charges a^ s 
and ay for Q 2 3> m 2 which is independent of the intermediate MS scheme: 
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where the coefficients are given by 
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It should be noted that this way of writing the two-loop coefficient r<iy in terms of ipy and 
ip v follows from the conformal ansatz. 

We now use the commensurate scale relation method to eliminate the scale ambiguity: 
Q is set to Q* using the single scale scale-setting approach |35| , such that all non-conformal 



terms proportional to ip v and ip v are absorbed into the running of the coupling |. This 
gives the next-to-leading order commensurate scale relation between a^ s and ay. To obtain 
the next-to-next-to leading order relation requires knowledge about the iV^-dependent part 
of the three-loop contribution. Thus we arrive at the following commensurate scale relation 
between a^ s and ay, 
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The next-to-leading order commensurate scale Q* is given by 
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which should be compared with the leading order commensurate scale Q* = 
y/sexp f — H + 2£ 3 J = 1.628-y/i. It should be noted that this way of writing the scale Q* 



differs slightly from the one used in |35[] in that it is written as the exponential of a partial 
fraction where the denominator is proportional to the \l/-function. This ensures that the 
scale Q* has sensible limits as ay — > or ay -^ oo (for Np = 3, Q* = 6.7Q in the limit 
ay — > oo). If the coupling freezes or is bounded for small scales, then the latter limit is of 
course not important. In addition Q* has the correct Abelian limit. The resulting commen- 
surate scale Q* is shown in Fig. |9] where it is also compared with the leading order scale. As 
can be seen from the figure, the next-to-leading order correction to the commensurate scale 
is small. The general convergence properties of the scale Q* as an expansion in ay is not 
known |36[ |. 



2 There also exists a multiple scale setting approach |T(| where one has different scales for each 
order of ay. However, for clarity we concentrate on only one of the procedures. In addition, as 
noticed in our earlier paper 
correct Abelian limit. 



13fl, the multiple scale setting procedure does not always have the 
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FIG. 9. The ratio of the commensurate scale Q* to y/s between the non-singlet width of the 
Z-boson and the heavy quark potential as a function of */s in next-to-leading (solid) and leading 
(dashed) order. 



The relation between a^iy/s) and ay{Q*) can be generalized to be valid for all scales, 
for which perturbation theory is applicable, by using the mass-dependent av(Q,rrii), 
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where the argument of a^ s is meant to indicated that the quark mass effects related to the 
running of the coupling are taken into account and m« being the pole-masses for the quarks 
which do not depend on Q. In addition we use the mass-dependent coupling ay(Q,rrii) 
and the mass-dependent coefficients of the ^y-function, ipy (Q,rrii) and ipy (Q, rrii), in the 
formula for Q* given by Eq. (|31~D . It should be noted that the scale Q* is only known to 
next-to- leading order. Similarly the evolution equation for ay(Q, rrii) is only known to next- 
to-leading order. Therefore we can only consistently use the next-to-leading order result 
when comparing with the treatment of finite quark mass effects in the MS scheme, i.e. 
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where the scale Q* should be the leading order result for consistency. 



C. Comparison 

Fig. (TO shows the relative difference between the next-to-leading order expressions for 
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NS 



in the MS and V schemes given by Eqs. (p6|) and (^) respectively. The predictions for 



(5), 



the width have been normalized to the same value at \/s = Mz using a^(Mz) = 0.118 and 
then evolved using next-to-leading order evolution in the respective schemes. 
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FIG. 10. The relative difference between the next-to-leading order expressions for o^ in the 
MS and V schemes respectively using next-to-leading order evolution. 



The comparison shown in Fig. |10| illustrates the relative difference between the predictions 
for a^ s in the MS and V schemes. In our earlier paper [13j we showed that the different 
ways of including the finite quark mass effects is smaller than ~ 0.1% by comparing the 
MS scheme with the analytic extension of the same which properly takes into account the 
flavor threshold effects analytically. Therefore the difference between the MS and V scheme 
predictions for a^ s can be attributed to the scheme dependence. This is illustrated by the 
fact that when using the next-to-leading order approximation for the commensurate scale, 
instead of the leading order one, the relative difference changes sign and even becomes larger. 
This sensitivity is a consequence of the scale dependence of the coupling, especially at small 
scales where the \l/-function is large. The proper inclusion of the finite quark mass effects is 
verified by the smoothness of the curve. 
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VI. SUMMARY AND CONCLUSIONS 

We have presented the calculation of ip v , the two- loop term in the Gell-Mann Low 
function for the ay scheme, with massive quarks. This gives for the first time a gauge 
invariant and renormalization scheme independent two-loop result for the effects of quark 
masses in the running of the coupling. Renormalization scheme independence is achieved by 
using the pole mass definition for the "light" quarks which contribute to the scale dependence 
of the static heavy quark potential. Thus the pole mass and the K-scheme are closely 
connected and have to be used in conjunction to give reasonable results. The results of 
the calculation are presented in numerical form due to the complexity of the integrations 
required. An important cross-check is the successful reproduction of the well-known QED 
results. 

The effective number of flavors in the two-loop coefficient of the Gell-Mann Low function 
in the ay scheme, Np V , becomes slightly negative for intermediate values of Q/m. This 
feature can be understood as anti-screening from the non-Abelian contributions and should 
be contrasted with the QED case where the effective number of flavors becomes larger than 
one for intermediate Q/m. For small Q/m the heavy quarks decouple explicitly as expected 
in a physical scheme, and for large Q/m the massless result is retained. 

The analyticity of the ay coupling can be utilized to obtain predictions for perturbatively 
calculable observables including the finite quark mass effects associated with the running 
of the coupling. By employing the commensurate scale relation method, observables which 
have been calculated in the MS scheme can be related to the analytic V-scheme without 
any scale ambiguity. The commensurate scale relations provides the relation between the 
physical scales of two effective charges where they pass through a common flavor threshold. 

As an example, we have shown how to calculate the finite quark mass corrections con- 
nected with the running of the coupling for the non-singlet hadronic width of the Z-boson 
compared with the standard treatment in the MS scheme. The analytic treatment in the 
V-scheme gives a simple and straightforward way of incorporating these effects for any ob- 
servable. This should be contrasted with the MS scheme where higher twist corrections due 
to finite quark mass threshold effects have to be calculated separately for each observable. 
The V-scheme is especially suitable for problems where the quark masses are important such 
as for threshold production of heavy quarks and the hadronic width of the r lepton. 
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APPENDIX A: ANALYTIC ow AT TWO LOOPS 



Finite quark masses are included naturally into the running of ay, thus providing an 
analytic definition of the gauge theory coupling. Furthermore, there is no scale ambiguity in 
a>v{Q) since the argument of the coupling is by definition the physical momentum transfer 
Q. These advantages can be carried over to the ordinary MS scheme by relating it to 
the physical ay scheme via a commensurate scale relation connecting the two schemes. It 
is in fact possible to combine the computational advantages of the MS scheme and the 
physical and analytic properties ay scheme into one common scheme, the analytic extension 
of the MS scheme [13|. However, as already mentioned in the introduction, the conformal 



coefficients in the commensurate scale relation between the ay and MS schemes does not 
preserve one of the defining criterion of the potential expressed in the bare charge, namely 
the non-occurrence of color factors corresponding to an iteration of the potential. This is 
probably an effect of the breaking of conformal invariance by the MS scheme. The breaking 
of conformal symmetry has also been observed when dimensional regularization is used as a 
factorization scheme in both exclusive [fT~5| , [T~6|1 and inclusive |17]] reactions. Thus, it does not 
turn out to be possible to extend the modified scheme 5ms beyond leading order without 
running into an intrinsic contradiction with conformal symmetry. For completeness we give 
the results of such an extension in this appendix. 



1. Commensurate scale relation between ay and a^ 



Our starting point for relating the ay and a^ schemes is the massless result for Q 2 3> m 2 



which is given by Eq. fl27|). Just as before we use the commensurate scale relation method 
to eliminate the scale ambiguity: the scale Q is set to Q* using the single scale scale-setting 
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between «ms and ay, 



This gives the following commensurate scale relation 
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The one- loop coefficient is the same as in our previous paper fT3f , but the two- loop one is 
changed due to the new result by Schroder [23||. However, the problem brought up in our 
previous paper regarding the anomalous contribution with a color factor proportional to 
CpNc is still there. This type of color factor corresponds to an iteration of the potential 
and thus cannot be part of the potential itself. The origin of this contribution is not clear, 
but it is probably an effect of the breaking of conformal invariance by the MS scheme. It 
should also be remarked that the conformal two-loop coefficient between the MS scheme and 
the ay scheme is large, indicating that there are large corrections between the two schemes. 
This is of great importance for observables like heavy quark production close to threshold 
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where the next-to-next-to leading order correction is known to be large in the MS scheme 
P7| . As another example the conformal coefficients for a^ s (s) in terms of a„g are also large, 
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where /i* is the commensurate scale between a^ s and the MS scheme. This should be 
compared with the conformal relation between a^ s and ay where the coefficients are not as 
large as indicated in Eq. (|30"D . Thus it is better to relate observables directly without using 
the intermediate analytic extension of the MS scheme. 

The commensurate scale Q* between the MS and V scheme is to next-to-leading order 
given by, 
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In the limit ay — > cxo the scale Q* becomes large (for Np = 3, Q* = 24Q in the limit 
ay — > oo). If the coupling freezes or is bounded for small scales, then this limit is of course 
not important. 

2. Definition of the Analytic a-^g 



The definition of the analytic a-^g is based on generalizing Eq. ( |Al| ) to be valid for all Q 
by using the mass- dependent ay(Q, rrii), 
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with rrii being the pole-masses for the quarks which do not depend on Q. In addition we 
use the mass-dependent coupling ay (Q, rrii) and the mass-dependent coefficients of the \Py- 
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In the above 



function, ip v '(Q, rrii) and i/j v (Q, rrii), in the formula for Q* given by Eq. 
definition we have only included terms to the order which we are working, i.e. next-to- 
leading order, since the effects from higher order terms on Q* is unknown. When going to 
even higher orders, the relation between the analytic a-^s and the ay scheme will contain 
large corrections as indicated in Eq. (|A1|) , reflecting the underlying large difference between 
the MS and ay schemes. 

We can also derive the \l/-function for «„g by taking the logarithmic derivative of Eq. (|A4|) 
with respect to Q. This gives 
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Re-expanding both sides of Eq. ( 
the first two terms in ^Ms(Q,m,, 



in ay using Eq. ( |A4|) and equating order by order gives 
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reflecting the non-trivial mass dependence of the ^-function. We note that when finite quark 
masses are included the first two terms in the ^-function are no longer universal but scheme 
dependent. 

We thus arrive at the following evolution equation for the analytic coupling a^^(Q,mi), 
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where ip^(Q,mi) an d "4 , j£g(Qi m i) are given by Eqs. (|A6|) and (|A7|), respectively, and rrii are 
the pole masses of the quarks. One complication which arises when solving the evolution 
equation is that the scale Q* has to be obtained recursively since Eq. ( |A3|) contains Q* also 
on the right hand side. In addition the approximation av(Q*,frii) = aMs(Q,mj) was used 
in the right hand side of Eq. ( |A3| ) when solving the evolution equation for aMs(Q, m i)- The 
evolution equation was solved for numerically using the classical Runge-Kutta algorithm. 
The resulting scale Q* calculated using Eq. ( |A3| ) is shown in Fig. |TT| . 
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FIG. 11. The ratio of the commensurate scale Q* to Q for the analytic extension of the MS 
scheme as a function of Q to leading (dashed) and next-to-leading (solid) order. 
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With the solution of the renormalization group equation for «ms we also obtain the \l/ 
function for the analytic extension of MS which is shown in Fig. 0. From the figure we also 



7 
6.5 

6 
5.5 

5 
4.5 

4 h 
3.5 



: : 






1 •. 




NT O 


\ 




LO 


- 


\ NLO 




- 


\> 






i 
i / 




LO step-function N F 


- 


T ^rrrr. 






1 


LO 


~~ *~ ^ — 


i — — _ 




- iim 


i 


. ; r ..,,, , mI 



10 



10 



10 10 

Q [GeV] 



FIG. 12. The normalized ^-function, —4>y[g(Q,mi)/(a^g(Q,mi)/Tr) in the analytic extension 
of the MS scheme compared to the conventional /3-function with discrete theta-function thresholds. 
The leading order results are shown as a dashed and dot-dashed curve respectively whereas the 
next-to-leading order results are shown as a solid curve and dotted curve respectively. 

see that the \l/-function in the analytic approach and in the massless step-function approach, 
with matching at the quark masses, follow each other closely except for small scales where 
they start to deviate. 



3. Comparing the Analytic a^ with a^g- 

We now compare the analytic 5ms with the conventional discrete theta-function treatment 
of flavor thresholds with matching at quark masses, a-^g. The relative difference between 



the two is shown in Fig. [13] both in leading and next-to leading order. 

As can be seen from the figure, the difference between the analytic and conventional 
treatment of quark masses in the running persists when going from leading to next-to-leading 
order. In fact we expect this difference to remain to all orders. The reason is that the \l/ 
function is not continuous in the massless approach and the stepsize at the thresholds is 
governed by the lowest order term ip(°\ Thus there will always be a finite difference between 
the continuous \l/-function and the one with theta-function thresholds. It is important to 
recognize that this feature is not eliminated by the fact that when going to higher orders 
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FIG. 13. The relative difference between the solutions to the 2-loop renormalization group 
equation using the analytic ^-function, &ms(Q), and conventional discrete theta-function thresh- 
olds, Ojjg(Q), with matching at quark masses. The solid curve shows the next-to- leading order 
result. For comparison the 1-loop result is shown as a dashed curve. The solutions have been 
obtained numerically starting from the world average [38] a-^(Mz) = 0.118. 



the dependence on the matching scale in the massless approach becomes smaller. The only 
way to include these mass effects in the ordinary MS treatment is by making a higher twist 
analysis to all orders in m 2 /Q 2 and Q 2 /m 2 . 
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